Coexisting charge modulation and ferromagnetism produces long 
period phases in manganites: new example of electronic soft 

matter. 



G.C. Milward, M.J. Calderon, P.B. Littlewood 

Cavendish Laboratory, Cambridge University, 
Madingley Road, Cambridge CBS OHE, UK 
(Dated: February 2, 2008) 



The phenomenon of colossal magnetoresistance in manganites is generally 
agreed to be a result of competition between crystal phases with different 
electronic, magnetic, and structural order; a competition which can be strong 
enough to cause phase separation between metallic ferromagnet and insulating 
charge modulated states . Nevertheless, closer inspection of phase 

diagrams in many manganites reveals complex phases where the two order 
parameters of magnetism and charge modulation unexpectedly coexist 0, 
Here we show that such experiments can be naturally explained within a 
phenomenological Ginzburg-Landau theory. In contrast to models where 
phase separation originates from disorder ja] or as a strain induced kinetic 
phenomenon [9], we argue that magnetic and charge modulation coexist in new 
thermodynamic phases. This leads to a rich diagram of equilibrium phases, 
qualitatively similar to those seen in experiment. The success of this model 
argues for a fundamental reinterpretation of the nature of charge modulation 
in these materials from a localised to a more extended "charge density wave" 
picture. The same symmetry considerations that favour textured coexistance 
of charge and magnetic order may apply to many electronic systems with 
competing phases. The resulting "Electronically soft" phases of matter with 
incommensurate, inhomogeneous and mixed order may be general phenomena 
in correlated systems. 



The manganese perovskites (RE^ X AE% + Mn0 3 , RE rare earth, AE alkaline earth) pro- 
vide a laboratory to study the interplay of a variety of magnetic, electronic and structural 
phases of matter in a strongly correlated electronic system. As in many strongly correlated 
electronic systems, the basic paradigm for manganite physics is the competition between the 
delocalising effects of the electron kinetic energy and the localising effects of the Coulomb 
repulsion, aided by coupling to lattice degrees of freedom. When the kinetic energy is dom- 
inant, one finds a metallic ground state with ferromagnetic alignment of the core moments. 
When the localising effects preponderate, instead we see charge and/or orbitally ordered 
ground states with substantial local lattice distortions from the near cubic symmetry of the 
metal, along with insulating behaviour and antiferromagnetism. One may tune between 
these two phases by many external parameters, especially chemical substitution, but also 
lattice strain, and magnetic field. The competition between metal and insulator is famously 
evident in the phenomenon of bulk colossal magnetoresistance, where a magnetic field tunes 
the conductivity of the material, and even more clearly in the strong tendency toward phase 
separation and inhomogeneity and regimes of percolative transport. 

The origin of charge and orbital ordered phases is still the subject of debate. Charge 
modulation has been traditionally seen as the ordering of Mn 4+ and Mn 3+ ions^J. More 
recently, the charge disproportionation of the Mn ions has been argued to be much smaller 
than one 1 131 ] but still the idea of two kinds of cation forming stripes prevails and 

is used to interpret the experiments. In such a scenario, one expects a density x of one 
kind of cation and 1 — x of the other. The charge modulation would be given by the 
averaged wave-vector q« (1 — a;)a* with a* the reciprocal lattice vector, aside from possible 
commensuration effects near special dopings (x = 1/2,2/3,3/4). 

However, this picture is not compatible with the experimental findings on commensu- 
rate and incommensurate modulation, summarised in FigO At half-doping, commensurate 
modulation is expected and found at low temperatures, together with antiferromagnetism of 
;he CE type |lo| . Above the Neel temperature, though, the modulation is incommensurate 



J HQ. 



When found to be charge modulated, underdoped samples (x < 0.5) do not 
show the relation q £a (l — x)a*, rather the modulation wave- vector is always commensurate 

samples (x > 0.5) show the ex pec ted incommensurate wave- vector below the Neel tem- 
perature, decreasing above it |20j]. However, no sign of discommensurations, but a uniform 



incommensurate modulation, has been found in recent experiments on Lai_j.Caa.MnO3 [21]. 

Another experimental conundrum is the coexistence of charge modulation and ferromag- 
netism despite their natural antipathy. For instance, at half-doping, the incommensurate 
modulation above the Neel temperature is accompanied by ferromagnetism [7]. A different 
electronic phase showing ferromagnetism and charge modulation has also been found at low 
temperature in Lan sCan sMn0 3 |6||. Slightly overdoped samples can also be ferromagnetic 



above Neel temperature 



22 



231 ] . Another example of coexistence is given by the underdoped 



0.3 < x < 0.5) Pr 1 _ a ,Ca a; Mn03. The ground state is commensurate and charge-modulated 
lfil \v\ but the antiferromagnetism is canted Q, 24 , |^ , showing a ferromagnetic compo- 
nent coexisting with the commensurate charge modulation. 

One might propose to explain these phenomena in terms of microscopic theory incorporat- 
ing the many different couplings and microscopic degrees of freedom, but this is a daunting 
task that might not be illuminating owing to its complexity. Here we propose a simple and 
more transparent phenomenological approach to this problem by means of Ginzburg-Landau 
theory. We will show how the close interplay of ferromagnetism and charge modulation con- 
spires to reproduce the experimental findings just discussed. 

Ginzburg-Landau theory allows the study of phase transitions in a phenomenological way 
and it consists in expressing the free energy as a power expansion of the order parameters 
and their gradients. The order parameters we consider here are the magnetisation M(r) 
and the charge-orbital modulation ip(r) = p(r)e i< - ( ^ c ' r+ ' ji< - r ^. r is the spatial coordinate, p is 
the amplitude of the modulation, Q c = — is a wave-vector commensurate with the lattice 
and (f> is the phase that would incorporate structures with incommensurate periodicities. To 
simplify the discussion we study a one dimensional scalar modulation, since charge modu- 
lation within a domain occurs only in one direction, n = 4 gives the correct periodicity for 
the lattice distortions measured in x = 0.5 as, though the charge modulation has period 2, 
the orbital order follows a zig-zag pattern with period 4. Notice that if V0 = 0, ip is a wave 
of amplitude p and wave-vector commensurate with the lattice. If V0 7^ 0, the wave-vector 
is Q c + (V0) and therefore, in general, cannot be expressed as a simple rational number of 
a*. 

The free energy density can be separated into three contributions: magnetisation, charge 



modulation and coupling terms. The first two are 




•F m = \a M {T - T C )M 2 + ^b M M 4 + ^(VM) 2 , (1) 

\a p (r - t co )p 2 + \b pP * + l -e P iy P ? + \e pP 2 cw - qo f + V cos^) (2) 



The magnetic energy, JF M , taken alone will describe a phase transition to homogenous 
magnetism below the Curie temperature Tc- F$ is the free energy extensively used to study 
commensurate- incommensurate phase transitions of charge density waves, spin density waves 
or modulated lattice distortions j^J. q a = 1/2 — x is the predicted deviation (by chemical 



a uniform incommensurate modulation with V0 = q a . On the other hand -r\p n cos(n0) is 
an Umklapp term that favours commensurability with <p = 2irj/n, j integer (for r] < 0). 
Taken alone, this describes two phases. Upon cooling below Tco, the amplitude p of the 
charge density wave rises from zero but provided n > 2, the Umklapp term is small and 
the modulation is incommensurate. As temperature is lowered, p grows, the Umklapp term 
may become dominant and a lock-in transition occurs if r\ is comparable to 

We now discuss coupling between the two order parameters. The lowest order coupling 
term which arises is d\p 2 M 2 so that there is a free energy penalty for homogeneous coex- 
istence of magnetism and charge modulation. Were this the only coupling term the free 
energy would be generally stabilised either by a homogenous magnetisation or by charge 
modulation, depending on which transition temperature is the larger. Next one can of 
course introduce uniform coupling terms of higher powers of M and p, but they make no 
qualitative changes unless they have a negative sign. More interesting is that there is a 
leading order coupling term in the gradient of the form d2p 2 M 2 {Wcf) — q Q ). The fact there 
is a term linear in the gradient is expected because there is no symmetry about x = 1/2; 
different signs of the gradient correspond physically to compression or extension of the CDW 
period, i.e. to extra "3+" or "4+" sites. One can also justify this term microscopically: if 
we consider the effect of charge modulation on the Fermi surface, then it is clear that if we 
choose a wave vector which does not match the chemical doping, one will be left with small 
pockets of carriers at the Fermi surface; these metallic electrons (or holes) are then available 
to mediate double exchange and thereby promote ferromagnetism. The asymmetry around 
x = 1/2 is due to the asymmetry between electron and hole pockets. Now note that this 
gradient term can be incorporated into Eq. |21by completing the square, and replacing q Q by 




The term 




favours 



qeff = qo -^= 1 -- X - d ^ . (3) 



The sign of c?2 is unknown a priori and we here choose it to be positive. Once this sign 
is fixed, however, this gradient coupling has profound consequences for the phase diagram. 
First, note that even if we are at commensurability (x = 1/2), if magnetism is present, then 
there is a tendency to incommensurate charge modulation. This reproduces the experiments 
of Chen et al. 0, on Lai/2Cai/2Mn03, where the onset of charge modulation is incommensu- 
rate, and accompanied by ferromagnetism - which is replaced by Neel order at the transition 
to the commensurate phase. The incommensurate phase of Pri^Cai^MnC^ (see Fig. [H) is 
paramagnetic but accompanied by the onset of ferromagnetic spin fluctuations 22 1- 
The second feature of this term is that if x < 1/2, it is possible for coexisting magnetism to 
"cancel" the chemical tendency to incommensurability, and we note that canted magnetism 



is generally reported |2 I 25] in the undordoped manganitos that show c-onmicnsurato 
charge modulation. No such cancellation is possible for x > 1/2, and thus the dog-leg 
dependence of wave- vector on doping shown in the inset to Fig. ^ is indicated. 

These features are reproduced by numerical minimisation of the coupled free energy with 
appropriate values of the parameters. Figure |2] shows the generic form of the phase diagram 
that can be obtained. In Fig. El we show an explicit evaluation of the temperature- and 
doping-dependence of the magnetism and (in) commensurability for parameters chosen to 
approximately reproduce the experimental regimes. The parameter range that can be used 
without altering the main features in Figs. El H3 is quite wide provided charge modulation 
dominates over magnetic order. 

Throughout much of the phase diagram shown in these figures, the order parameters M, 
p and V0 are approximately uniform in space, at least for the parameters we have cho- 
sen. However, close to the commensurate-incommensurate transition, the phase modulation 
becomes non uniform, and the phase gradient is built up by periodic discommensurations 
where the phase advances through 2n/n. In such an inhomogeneous state, magnetism is nat- 
urally enhanced at the boundary and the amplitude of the charge modulation suppressed 
(see Fig. H}. Another situation where inhomogeneity is enforced is at a magnetic domain 
wall, where it can be energetically preferable to have a sharp wall stabilised by an insertion 
of local charge modulation. Such a phenomenon might be responsible for the anomalously 



large resistance reported for magnetic domain walls in LaCaMnOa [281 129]. 

When interpreting experimental results in the light of the theory presented here, a couple 
of issues must be kept in mind. Firstly, the theory proposed is a mean field theory that 
cannot describe strong fluctuations as its solutions are uniformly ordered or disordered 
phases. However, experimentally there are some regimes, like the incommensurat e p hase 
of PTi/ 2 Ga 1 / 2 Mii0 3 , where strong ferromagnetic fluctuations have been measured [la, \27 \. 
rather than the long range magnetic order we would predict. Secondly, in order to keep the 
calculations tractable, the phase transitions have been forced to be continuous. This explains 
why the reentrant magnetism above Tl shown in Fig. El appears only when Tq > Tco- 
If this theory were generalised to include discontinuous phase transitions, this condition 
would relax. Another consequence of assuming continuous phase transitions is that phase 
separation cannot be predicted. In real systems phase separation is possible since strain 
a, ly, or disorder ja], can make more or less localised phases dominate within a given region. 
Orbital ordering is described by a vector order parameter, thus our simple model cannot 
address the complexity of different orbitally ordered phases that have been proposed jl^ . 

The Ginzburg-Landau phenomenology we propose is capable of systematizing some puz- 
zling data for manganites near x = 1/2, but of course the propensity for mixed and homoge- 
neous phases is driven by the underlying physical parameters that make the energetic cost 
of spatial fluctuations low. This "electronic softness" means that as well as spatially disor- 
dered "phase separation" , one finds new ordered phases which are long period arrangements 
of the two competing orders. It may indeed be that this potential for textured electronic 
phases is a hallmark of electronic oxides near the Mott transition J3C| . seen perhaps in the 
co-existence of density waves and superconductivity in the cuprates [3l|. 

Acknowledgements. 

The authors thank L. Brey and N. Mathur for helpful discussions. PBL thanks the National 
High Magnetic Field Laboratory of the Los Alamos National Laboratory for hospitality. 
MJC acknowledges Churchill College, University of Cambridge, for the award of a JRF. 
This work was supported by the EPSRC and through the EPSRC Magnetic Oxide Network. 



[1] M.B. Salamon and M. Jamie, The physics of manganites: structure and transport, Reviews 
of Modern Physics 76 583 (2001) 



[2] N.D. Mathur and P.B. Littlewood, Mesoscopic textures in manganites, Physics Today 56, 
25-30 (2003). 

[3] N.D. Mathur and P.B. Littlewood, The self-organised phases of manganites, Solid State Com- 
munications 119, 271-280 (2001). 

[4] E. Dagotto, Nanoscale phase separation and colossal magnetoresistance, Springer Series in 
Solid State Sciences, Vol 136. 

[5] M. Uehara, S. Mori, C.H. Chen and S.-W. Cheong, Percolative Phase Separation underlies 
colossal magnetoresistance in mixed- valence manganites, Nature 399, 560-563 (1999). 

[6] J.C. Loudon, N.D. Mathur and P.A. Midgley, Charge-ordered ferromagnetic phase in 
La .5Cao.5Mn0 3 , Nature 420, 797-800 (2002). 

[7] C.H. Chen and S.-W. Cheong, Commensurate to Incommensurate Charge Ordering and Its 
Real-Space Images in Lao.5Cao.5Mn03, Physical Review Letters 76, 4042 (1996). 

[8] A. Moreo, M. Mayr, A. Feiguin, S. Yunoki, and E. Dagotto, Giant Cluster Coexistence in 
Doped Manganites and Other Compounds, Physical Review Letters 84, 5568-5571 (2000). 

[9] K.H. Ahn, T. Lookman and A.R. Bishop, Strain-induced metal-insulator phase coexistence in 
perovskite manganites, Nature 428, 401-404 (2004). 
[10] J. Goodenough, Theory of the Role of Covalence in the Perovskite- Type Manganites [La, 

M(H)]Mn03, Physical Review 100, 564 (1955). 
[11] J. Garcia, M.C. Sanchez, J. Blasco, G. Subfas and G. Proietti, Analysis of the x-ray resonant 
scattering at the Mn K edge in half-doped mixed valence manganites, Journal of Physics: 
Condensed Matter 13, 3243-3256 (2001). 
[12] L. Brey, Continuous Charge Modulated Diagonal Phase in Manganites, Physical Review Let- 
ters 92 127202 (2004). 
[13] Coey M, Charge ordering in oxides, Nature 430 155-156 (2004). 

[14] S. Mori, T. Katsufuji, N. Yamamoto, C.H. Chen and S.-W. Cheong, Microstructure related 

to charge and orbital ordering in Pro.5Cao.5Mn03, Physical Review B 59, 13573 (1999). 
[15] R. Kajimoto, H. Yoshizawa, Y. Tomioka and Y. Tokura, Commensurate-incommensurate 

transition in the melting process of orbital ordering in Pro.sCao.sMnOs: A neutron diffraction 

study, Physical Review B 63, 212407 (2001). 
[16] M.v. Zimmermann, C.S. Nelson, J.P. Hill, D. Gibbs, M. Blume, D. Casa, B. Keimer, Y. 

Murakami, C.-C. Kao, C. Venkataraman, T. Gog, Y. Tomioka, and Y. Tokura, X-ray resonant 



scattering studies of orbital and charge ordering in Pri-^Ca^MnOs, Physical Review B 64, 
195133 (2001). 

[17] S. Larochelle, A. Mehta, N. Kaneko, P.K. Mang, A.F. Panchula, L. Zhou, J. Arthur, and M. 
Greven, Nature of e g electron order in Lai_ x Sri + - E Mn04, Physical Review Letters 87, 095502 
(2001). 

[18] Z. Jirak, S. Krupicka,Z. Simsa, M. Dlouha and S. Vratislav, Neutron diffraction study 
of Pri-zCa^MnOs perovskites, Journal of Magnetism and Magnetic materials 53, 153-166 
(1985). 

[19] C.H. Chen, S.-W. Cheong, and H.Y. Hwang, Charge ordered stripes in Lai-^Ca^MnOs, Jour- 
nal of Applied Physics 81, 4326 (1997). 

[20] C.H. Chen, S. Mori and S.-W. Cheong, Anomalous melting transition of the charge ordered 
state in manganites, Physical Review Letters 83, 4792-4795 (1999). 

[21] J.C. Loudon, S. Cox, A.J. Williams, J.P Attfield, P.B. Littlewood, PA. Midgley and 
N.D. Mathur, Uniform periodicity in stripe phase Lai-^Ca^MnOs at arbitrary doping, 
cond-mat/0308581 , submitted to Physical Review Letters. 

[22] P. Schiffer, A. P. Ramiresz, W. Bao and S-W. Cheong, Low Temperature Magnetoresistance 
and the Magnetic Phase Diagram of Lai_a;Ca x Mn03, Physical Review Letters 75, 3336 (1995) 

[23] R. Kajimoto, H. Yoshizawa, H. Kawano, H. Kuwahara, Y. Tokura, K. Ohoyama and M. 
Ohashi, Hole-concentration-induced transformation of the magnetic and orbital structures in 
Ndi-^Sr^MnOa, Physical Review B 60, 9506 (1999). 

[24] Y. Tomioka, A. Asamitsu, H. Kuwahara, Y. Moritomo, and Y. Tokura, Magnetic-field-induced 
metal-insulator phenomena in Pri-^Ca^MnOs with controlled charge-ordering, Physical Re- 
view B 53, 1689 (1996). 

[25] Yoshizawa H, Kawano H, Tomioka Y, Tokura Y, Neutron diffraction study of the magnetic- 
field-induced metal-insulator transition in Pro.yCao.3Mn03, Physical Review B 52, R13145- 
R13148 (1995). 

[26] J.C. Toledano and P. Toledano, The Landau Theory of phase transitions. World Scientific, 
1987 

[27] R. Kajimoto, T. Kakeshita, Y. Oohara, H. Yoshizawa, Y. Tomioka, and Y. Tokura, Anomalous 
ferromagnetic spin fluctuations in an antiferromagnetic insulator Pri_ x Ca x Mn03, Physical 
Review B 58, R11837 (1998). 



[28] Mathur ND, Littlewood PB, Todd NK, Isaac SP, Teo BS, Kang DJ, Tarte EJ, Barber ZH, 
Evetts JE, Blamire MG, Resistance of a domain wall in Lao.7Cao.3MnO"3, Journal of Applied 
Physics 86, 6287-6290 (1999). 

[29] Rzchowski MS and Joynt R, Electronic inhomogeneity at magnetic domain walls in strongly- 



correlated systems, cond-mat/0304287 
[30] S.A.Kivelson E. Fradkin and V.J.Emery, Electronic liquid-crystal phases of a doped Mott 

insulator, Nature 393, 550-553 (1998). 
[31] K. McElroy et al, Homogenous nodal superconductivity coexisting with inhomogeneous charge 



order in strongly underdoped Bi-2212, cond- mat/0404005 



u 

O 
-»-» 
o 

> 

> 
c3 



0.35 - 



0.2 



i i i 




i i i 



Pr 0.6 Ca 0.4 



Pr f^a 

ri l/2^ d l/2 



La 0.45 Ca 0.5: 



La 3/8 Ca 5/8 



Pr 

ri l/3^ d 2/3 



0.5 

0.3 ho.4 1- 



o.3 - q/a*~l-x* 



0.3 0.5 0.7 

■ ■ I X ■ ■ ■ 




La 1/3 Ca 2/3 



20 150 200 250 

Temperature(K) 



290 



FIG. 1: Wave vector of the modulation q/a* versus temperature for Pr i T C& x MnO.i and 



Lai_ x Ca x Mn03 at different dopings (x > 0.5 data taken from Figs. 1 and 2 in 20 and x = 0.4 
taken from 3]). The same kind of behaviour has been reported for LaipCai^MnO^ 2] (not shown 
here). At low temperatures q/a* ~ (1 — x), as shown in the inset, and decreases above the Neel 



□ 



191 ] . The circles correspond to the low 



temperature Taf- The inset shows q/a* versus x as in 
temperature values of the curves in the main panel while the squares are data for Lai-^Ca^MnOs 
taken from 21]. 
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FIG. 2: Schematic phase diagram which results from the minimisation of the free energy. The scale 
of the axes depends on the particular parameters used. The commensurate order phase just above 
x = 0.5 has not been observed but we predict it can be relevant for dopings very close to x = 0.5 
for highly insulating manganites. The complex phases arise provided Tc > Tco, a condition 
that can be relaxed if the model is extended to account for discontinuous phase transitions. The 
values of Tqo an d Tc are direct parameters in the model whereas the lock-in temperature Tl is a 
consequence of the competition between the Umklapp and incommensurate modulation terms. 








FIG. 3: Results within Phase Modulation Approximation for a particular choice of parameters. 
These complement and clarify the schematic view in Fig. |2I Top: magnetisation in the doping- 
temperature plane. The low temperature phases are non magnetic for x > 0.5 (corresponding to 
antiferromagnetism in experiments) and weakly ferro mag netic for x < 0.5 (corresponding to the 
canted antiferromagnetism found in Pri_ x Ca x Mn0 3 QflQ)- 

Above the Neel temperature 



there is a reentrant magnetisation that has been reported for x > 0.5. Bottom: deviation from 
commensurability 1/2— q in the doping-temperature plane. The surface consists of three planes that 
correspond to commensurate order V</> = for the ground state of the charge ordered underdoped 
samples, incommensurate order for the overdoped samples that follow 1/2 — x as shown in the inset 
of Fig. ^ and a change of wave vector with temperature above the Neel transition as shown in the 
main panel of Fig^ 
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FIG. 4: Example solutions of the total Free Energy when all the order parameters are allowed to 
change spatially, (a) The magnetisation (green) decays in the centre of a magnetic domain wall 
leading to the appearance of charge modulation (red), (b) In a discommensuration, the phase of 
the charge modulation (blue) changes by vr/2, the amplitude of the charge modulation (red) is 
suppressed, and the magnetisation (green) enhanced. 



